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ABSTRACT 

Galaxy counts and recent measurements of the luminosity density in the near-infrared (NIR) have 
indicated the possibility that the local universe may be under-dense on scales of several hundred 
megaparsecs. The presence of a large-scale under-density in the local universe could introduce signif- 
icant biases into the interpretation of cosmological observables, and, in particular, into the inferred 
effects of dark energy on the expansion rate. Here we measure the K—hanA luminosity density as a 
function of redshift to test for such a local under-density. We select galaxies from the UKIRT Infrared 
Deep Sky Large Area Survey (UKIDSS-LAS) and use spectroscopy from the Sloan Digital Sky Survey 
(SDSS) and other redshift surveys to generate a i^— selected catalog of ~ 35, 000 galaxies that is 
> 90% spectroscopically complete to K^s, = 16.3. The relative depth and wide area of our sample 
allow for a detailed measurement of the iiT— band luminosity density as a function of redshift in the 
range 0.005 < z < 0.2. We find that the overall shape of the z = rest-frame K—h&nA luminosity 
function [M* = —21.38 ± 0.04 and a = —1.02 ± 0.03) appears to be relatively constant as a function 
of environment and redshift out to z ^ 0.2. We find a local {z < 0.07) luminosity density that is 
in good agreement with previous studies. At z > 0.07 we detect a rising luminosity density, and at 
z > 0.1, it is roughly ~ 1.5 times higher than that measured locally. This suggests that the stellar 
mass density as a function of redshift follows a similar trend. Assuming that luminous matter traces 
the underlying dark matter distribution, this implies that the local mass density of the universe may 
be lower than the global value on a scale and amplitude sufficient to introduce significant biases into 
the determination of basic cosmological observables, such as the expansion rate. An under-density of 
roughly this scale and amplitude would be sufficient to resolve the apparent tension between direct 
measurements of the Hubble constant and those inferred by Planck. 
Subject headings: cosmology: observations — galaxies: fundamental parameters 



1. INTRODUCTION 

The universe is generally assumed to be isotropic 
and homogeneous on very large scales. This allows 
for the development of cosmological models and obser- 
vations to constrain those models, provided that the 
observations are made over a sufficiently large volume 
to average over so called "cosmic variance", or sys- 
tematic measurement biases due to large-scale struc- 
ture. Observations of the cosmic microwave back- 
ground (CMB) indicate the universe was very homo- 
geneous at z ~ 1100 and measurements of the kine- 
matic Sunyaev-Zeldovich (kSZ) effect appear to indicate 
that the present day universe is homogeneous on scales 
greater than ^ 1 Gpc (iGarcia-Bellido fc Haugb0llell2OO8l : 
iZhang fc Stebbinsll20lH ). However, the large-scale ho- 
mogeneity of the universe on smaller scales has not been 
measured directly. 

Observed large-sc ale structures, suc h as the > 400 Mpc 
Sloan Great Wall (jGott et al.ll2005[ ). demonstrate the 
existence of inhomogeneity on very large scales. In 
addition, it has been shown that voids on similarly 
large scales may explain the "cold spots" in the CMB 
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(llnoue fc SilU 120061 ) ■ While such structures had previ- 
ously been thought to be departures from the e xpecta- 
tion of ACDM cosmology (|Sheth fc Diaferioll201ll ). it has 
now been found that such large-scale structures are not 
only to be expected in the current concordance cosmol- 
ogy, but also that even larger structur es are likely to be 
discovered as survey volumes increase (jPark et al.ll2012f) . 
Recent cosmological modeling efforts have demon- 
strated that large-scale structure in the local universe 
may introduce significant systematic errors into the 
measurement of cosmological observables and hence 
the interpretation of these observables in the context 
of a given cosmolo g ical model (for recent reviews see 
iBoleiko et all 120111: iClarksonI 120121 ). In particular, 
so called "void models" place the observer inside 
a large local under-density to provide for the ap- 
parent acceleration of the expansion of the universe 
Moffa t fc Tatarski 1992, 199 5; Cc l cricr 2000; T omite 
llguchi et all 120021: lAlnes et all 
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Chung fc Romano! 120061: Enqvist fc MattssonI 
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12007 
Yoo et all 120081: iGarcia-B dhdo fc Haugb0llc 200^" 
'Alexa nder et al.ll2009HGarcia-Bellido fc HauaboUc 200", 
Febru ary et al.l 120101: iCelerier et aL feolO; Biswa s et al.| 
12010; Marra fc Paakkoncn 2010; Clarks on fc MaartensI 
12010; Bolc iko fc Sussman 2011: Nishikawa et al.l 120121: 
Valkcnburglloi^! 

The basic idea underlying these void models is that if 
an observer lives near the center of a large under-density, 
then that observer will witness a local expansion of the 



universe that is faster than the global expansion. This 
would result in a locally measured Hubble constant that 
is larger than the global expansion rate and appear ob- 
servationally as an accelerating expansion. 

In their current form, void models without a cos- 
mological constant do not appear to be viable alter- 
natives to dark-energy-dominated universes, as they 
have problems in simultaneously fitting all cosmo- 
logical obs ervable s ([G arcia-Bellido & Haugb0llc 2008; 
Zibin et a ll 120081: iMoss ct al. 2011; Zhang & Stebbina 
2011; Ri ess et al.l 2011; Zumalacarrcgui et al. 2012). 
However, the exploration of these types of cosmologi- 
cal models — and other models which explore the effects 
of large-scale inhomogeneity — has highlighted the need 
for a more thorough understanding of extremely large- 
scale structure in t he lo cal univ erse JM arra & Notarij 
20ll iMar ra et al." 2012, 201 1 iBull fc Clifton 2012!; 



Mishra et al. 2012 . 2013: Valkenburg et al. 2012 . 2013) . 

In particular , "m i nimal void" scenarios (e.g., 
[Alexander et al.l 120091: iBoleiko fc SussmanI I2011D have 
shown that very simple models placing the observer 
near the center of an under-density that is ~ 250 h~^ 
Mpc in radius (extending to z ^ 0.08) and roughly half 
the density of its surroundings are sufficient to explain 
the acceleration observed via type la supernovae. While 
these models are simplistic, they make clear that an 
observer's location with respect to large-scale structures 
could have profound implications for that observer's 
measurement of cosmological observables. 

In this study, we wish to make an estimate of the mass 
density of the nearby universe as a function of redshift 
to test for large-scale inhomogeneity, and in particular, 
for a large local under-density. While we cannot directly 
probe the underlying dark matter distribution, we can 
make a robust estimate of the near-infrared (NIR) lumi- 
nosity density, whic h is a good tracer of the overall stellar 



mass dens ity (e. g.. [deJondll996HBell fc de Jon3 120011: 
iBeU et al.l [201031 kirbv et al.ll2008[ ). The stellar masses 



of galaxies, on average, have been shown to be corre- 
lated with the ma ss of their host dark matter halos (e.g., 
iWang et al.ll2012l ). Furthermore, simulations have shown 
that on much larger (~ 100 Mpc) scales at low redshift, 
the spatial distribution of baryons should be an unbi- 
ased tracer of the underlying dark matter distribution 
(jAngulo et al.l [20131 ). 

Thus, in terms of the linear bias parameter, here we 
will assume 6 = 1 in the relation Sg = bS, where Sg is 
the relative density contrast given by the distribution of 
galaxies, i.e., {pstars.z - P star s)Ip stars ^ ^nd d is the total 
mass density contrast, {pz—'p)/'p. At optical wavelengths, 
the linear bias parameter h as been ob served to approach 
unity at low redshifts f e.g.. | Marinoni e t al. 2005) , and in 
the NIR, Mailer et al.l ()2005| ) have measured a value in 
the X— band of hxiz « 0) = 1.1 ± 0.2, suggesting the 
NIR bias parameter follows a similar trend. 

Therefore, a measurement of the NIR luminosity den- 
sity can provide an estimate of the underlying mass den- 
sity of the universe. Likewise, a measured change in 
NIR luminosity density as a function of redshift could 
signal a corresponding change in the underlying total 
mass density. Thus, in the NIR at low redshifts, where 
dust extinction is minimal and iC— corrections are small 
and nearly independent of galaxy type, statistical studies 
of galaxies provide a means of probing local large-scale 



structure. 

Several studies of NIR galaxy counts have found that 
the local space density of galaxies appears to be low 
by '-^ 25 — 50% compared to the density at distances 
of ~ 300 h~^ Mp c or b e yond z ^ 0.1 (iHuang et al" 



01 ^ dUU n " Mpc or be yong z ^ U.l (IMuang et a 
19971: iFrith et all IMil (20051 : IBusswell et al.l 1200 



Keenan et al.l l2010al) ^ If the space density of galaxies 



is rising as a function of redshift, then a corresponding 
rise in NIR lumino s ity de nsity should also be present. 

In iKeenan et al.l (|2012l ). we probed the NIR luminos- 
ity function (LP) just beyond the local volume at red- 
shifts of z ^ 0.2. We found that the NIR luminosity 
density at z ^ 0.2 appears to be ~ 30% higher than 
that measured at 2 ~ 0.05. This measured excess could 
be considered a conservative underestimate, because we 
avoided known over-densities, such as galaxy clusters, in 
our study. We note that our result cannot be consid- 
ered conclusive, given possible systematics due to cos- 
mic variance in our measurement. However, taken in the 
context of other measurements of the NIR luminosity 
density from the literature, our result is consistent with 
other studies that found a higher luminosity density at 
z > 0.1 than that found locally. Furthermore, we showed 
that all measurements of the luminosity density could be 
considered roug hly consistent with the voi d radial den- 
sity profile from IBoleiko fc SussmanI ()201lD . which they 
claim could mimic the apparent acceleration of the ex- 
pansion of the universe. 

In the past, wide-area NIR surveys have generally not 
gone deep enough to probe beyond the very local uni- 
verse, and deep surveys have been carried out over rela- 
tively small solid angles, such that good counting statis- 
tics cannot be achieved at lower redshifts. The result is 
that, in a study such as that performed in IKeenan et al.l 
(2012), we are comparing measurements made with quite 
different photometry and methodology, and the results 
may suffer from unknown biases. 

Recently, however, the completion of the relatively 
deep and wide UKIR T Infrared Deep Sky Large Area 
Survey fUKIDSS-LAS. iLawrence et al.l[2007[ ). combined 
with s pectroscopy from the Sloan Digital Sky Survey 
(SDSS, l^rFetSlIlQQO.), the Two-degree Field G alaxy 
Redshift Survey (2dFGRS, ICoUess et al.) l200l . and 
other spectroscopic data from the archives, provides for 
a NIR-selected spectroscopic sample of galaxies that is 
both wide enough on the sky and deep enough pho- 
tometrically to sample a relatively broad range in red- 
shifts in the nearby universe with good statistics. In 
this paper, we draw a i^— band selection of galaxies from 
the UKIDSS-LAS that is highly spectroscopically com- 
plete. We then study the ii'— band galaxy LF over a 
large volume in the redshift range 0.005 < z < 0.2. With 
these data we are able, for the first time, to measure 
the X— band galaxy LF as a function of redshift using 
consistent photometry and methodology. 

The structure of this paper is the following. We dis- 
cuss the sample selection in Section [2j We estimate the 
i^— band luminosity density as a function of redshift and 
possible biases in our measurements in Section [3l We 
summarize our results in Section |4l Unless otherwise 
noted, all magnitudes listed here are in the AB magni- 
tude system (toab = 23.9 — 2.5 log^o ff with /^ in units 
of /iJy). We assume Hm — 0.3, ft\ = 0.7, and h = 0.7 in 
our calculations. 
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Figure 1. The region on tlie slcy wliere tlie UKIDSS-LAS DR8 and tiie SDSS DR9 overlap and tiie spectroscopic completeness is > 90% 
(~ 536.9 deg^ total). Tlie black rectangles represent 1° X 1° areas of high completeness within the overlap region. The boxes drawrn in 
green, blue, and orange denote the three subregions in this sample that arc subsequently dealt writh separately to address possible biases 
and cosmic variance in this study. 



2. CATALOG GENERATION 

To create a stellar mass selected sample, we have re- 
trieved data from the WFCam Science Archive (WSA) 
for all galaxies with NIR photometry in the UKIDSS- 
LAS DR8 in the iC-band (Kyega < 14.4, which is equiv- 
alent to Kab < 16.3), where overlap with the SDSS, the 
2DFGRS, and other surveys provides for highly com- 
plete spectroscopy of bright galaxies (> 90% complete 
for Kab = 16.3). This region of sky covers 536.9 deg^, 
and the selection contains 35, 342 galaxies at a median 
redshift of z ~ 0.1. 

We then used the SDSS-III Sky Server CasJobsQ inter- 
face to retrieve optical photometry and redshifts from 
the SDSS DR9 over the areas covered in all bands 
(F, J, H, and K) in the UKIDSS-LAS. We then cross 
correlated the positions of the if —selected objects with 
those from the SDSS using a search radius of 2". 

A small fraction (< 1%) of objects identified in the 
UKIDSS catalogs did not have a counterpart within < 
2" in the SDSS. We found that these objects generally 
appeared to be spurious iiT— band detections (checking by 
eye in the imaging data from UKIDSS), so we excluded 
these objects from the final catalog. Thus, all of the 
UKIDSS X— band selected objects in our catalog have 
counterparts classified as primary target galaxies in the 
SDSS. We find that 95% of UKIDSS objects matched 
to SDSS counterparts have angular separations of < 0'.'5, 
and the average separation between UKIDSS objects and 
their SDSS counterparts is 0'.'2. 

We also downloaded the 2dFGRS redshift catalogs 
from the VizieR online service^. We cross-matched 2dF 
objects to their UKIDSS counterparts within a radius of 
2". Cross-correlation with the 2dFGRS significantly im- 
proved overall completeness in the 2dF equatorial region 
from lO'' < RA < 15'' (subregion 2 denoted in Figure[T]). 
We also cross-correlated with published redshifts in the 
NASA-Sloan Atlaffl and the NASA/IPAC Extragalactic 
Database^ to augment the completeness of the sample. 

^ skyserver.sdss3.org/CasJobs 

® http://vizier.cfa.harvard.edu/viz-bin/VizieR-3 

7 www.nsatlas.org 

® ned.ipac.caltech.edu 
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Figure 2. Completeness as a function of /C— band apparent mag- 
nitude. The average completeness of the full sample of 35, 342 
galaxies with K\b < 16.3 is ~ 95%. 

We restricted our final catalog to regions on the sky 
which are > 90% spectroscopically complete at Kab < 
16.3 in order to minimize possible biases associated with 
the fact that our sample is selected in the K—h&nd, while 
the targets for the surveys providing the redshifts for this 
work were primarily optically selected. This resulted in 
a catalog of 35, 342 galaxies over the 536.9 deg^ region 
shown in Figurc[Tl after excluding stars from the catalog, 
as we describe in the following subsection. We show the 
completeness of the sample as a function of apparent 
magnitude in Figure [5] and the redshift distribution in 
Figure [31 

2.L Star- Galaxy Separation 

We investigated the g— J vs. J—K color-color diagram 
for objects in our catalog to determine the reliability of 
the star-galaxy classifiers offered in the WSA and SDSS 
archives (UKIDSS "mergedclass" and SDSS "type")- We 
show a gJK color-color diagram for 3" aperture mag- 
nitude colors in Figure [H Orange points show objects 
classified by morphology as extended in both UKIDSS 
and the SDSS (UKIDSS mergedclass = -1 and SDSS 
type = 6). Light blue points show objects classified as 
point sources in both UKIDSS and the SDSS (UKIDSS 
mergedclass = 1 and SDSS type = 3). Red points show 
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Figure 3. The redshift histogram for the spectroscopic sample of 
galaxies in this study. Subrcgion 2 overlaps with the Sloan Great 
Wall identified bv iGott ct al.. (,2005l ) . and the peak in the redshift 
histogram in this region near z ~ 0.08 is primarily due to this 
structure. 
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Figure 4. A 3" aperture magnitude color-color separation for 
all objects at ii" < 16.3 in the raw sample drawn from the 
WSA and SDSS Skyserver. Orange points show objects classi- 
fied by morphology as extended in both UKIDSS and the SDSS 
(UKIDSS mergeddass = -1 and SDSS type = 6). Light blue 
points show objects classified as point sources in both UKIDSS 
and the SDSS (UKIDSS mergeddass = 1 and SDSS type = 3). 
Red points show spectroscopically identified galaxies. Purple stars 
show spectroscopically identified stars. The dashed line shows 
the best color separation boundary between stars and galaxies: 
g - J = 5.28 X [J - K] + 2.78. Approximately ~ 2% of extended 
sources lie above the separation boundary and ~ 0.5% of po int 
sources lie below. Based on analyses described in Section 12.11 we 
conclude that by excluding all sources above the separation bound- 
ary, as well as all sources identified by morphology as point sources 
in the SDSS and WSA (regardless of color), we can achieve a star- 
galaxy separation that is robust at better than the 1% level. 



spectroscopically identified galaxies. Purple stars show 
spectroscopically identified stars. The dashed line shows 
the best color separation boundary between stars and 
galaxies: g — J = 5.28 x [J — K] + 2.78. Approximately 
~ 2% of extended sources lie above the separation bound- 
ary and ~ 0.5% of point sources lie below. 

Next, we looked at the images in the SDSS Sky Server 
database of the point sources that lie below the separa- 
tion boundary and the extended sources that lie above. 
For these point sources, it was clear from the SDSS imag- 



ing that they are point-like in morphology. We expect 
that most of these objects are stars with rare colors or 
bad photometry in one band (with the exception of a few 
quasars, which are not of interest for this study). 

The vast majority of the objects classified as extended 
that lie above the separation boundary lie along the stel- 
lar main sequence for this color separation. Upon in- 
vestigating the imaging for these sources, we found that, 
while there are definitely a few galaxies among them, the 
vast majority (~ 99%) look like point sources that are 
either smeared out a bit or blended with another object. 
Thus, we conclude that the majority of objects that are 
classified as extended but lie above the boundary are, in 
fact, stars. 

We note that 15 spectroscopically confirmed galaxies 
lie above the separation boundary. Upon closer inspec- 
tion of SDSS imaging for these sources, we found that 
roughly a third appeared to be point-sources on the stel- 
lar main sequence (likely stars misclassified as galaxies 
due to noisy spectra). Another third appeared to be 
galaxies blended with a star, and the final third appeared 
to be galaxies by morphology, but with stellar colors. 
There are only a handful of spectroscopically confirmed 
stars that lie below the separation boundary, and all ap- 
pear to be point sources. 

Based on these analyses, we conclude that by excluding 
all sources above the separation boundary described in 
Figure |4l as well as all sources identified by morphology 
as point sources in the SDSS and WSA (regardless of 
color), we can achieve a star-galaxy separation that is 
robust at better than the 1% level. 

2.2. Petrosian Aperture Clipping in UKIDSS 
Photometry 

The sky subtraction algorithm in the pipeline for 
UKIDSS photometry is such that there exists an upper 
limit of 6" on the Petrosian aperture radius (correspond- 
ing to a circular aperture with a radius of 12"). This 
causes the total flux from galaxies that subtend large 
solid angles to be underestimated. This implies that if 
we use Petrosian magnitudes, then for some galaxies we 
will underestimate their luminosities, while other galax- 
ies will be lost entirely from the sample, because the 
aperture clipping has pushed them to a fainter apparent 
magnitude than the selection limit of the sample. 

We retrieved the i^— band Petrosian aperture radii for 
our sample from the WSA and found that roughly ^ 10% 
of galaxies in the sample had their Petrosian apertures 
clipped at 6". We show the fraction of galaxies for which 
the Petrosian aperture was clipped at 6" as a function of 
redshift in Figure [5] Clearly, the underestimation of total 
flux is a much stronger effect at low redshift (z < 0.1). 
Thus, in general, brighter (larger) galaxies at low redshift 
are most strongly affected by Petrosian aperture clipping. 

Instead of omitting galax ies for which the Petrosian 
apertures were clipped (^e.g.. lSmith et al.ll2009() . we have 
devised a method to compensate for the effects of the un- 
derestimation of flux. Short of redoing all the photom- 
etry, there is no way to know exactly what fraction of 
light was lost due to the aperture clipping. However, for 
each galaxy the WSA provides photometry for circular 
apertures ranging in size from 1" to 12" in radius. Using 
these measurements, we determined a rough light curve 
for each galaxy affected by Petrosian aperture clipping. 
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Figure 5. The fraction of galaxies with their Petrosian aperture 
clipped at 6" as a function of redshift in the UKIDSS-LAS. We 
correct for this effect at z < 0.1 by extrapolating the light curve 
derived fr om a range of circular aperture magnitudes, as described 
in Section [22] 

We then extrapolated this light curve out to the SDSS 
z— band Petrosian radius provided for each galaxy. In the 
NASA-Sloan Atlas (NSA), new and improved photome- 
try is provided for SDSS galaxies at z < 0.055. In the 
case of galaxies in our sample that are included in the 
NSA catalog, we extrapolate to the new z— band Pet- 
rosian radius provided by the NSA. 

At z > 0.1 the aperture clipping problem is minimal 
(see Figure [5]). Given the nature of this study to investi- 
gate the possibility of a local (z < 0.1) under-density, we 
elected only to correct for aperture clipping at z < 0.1, 
to be as conservative as possible in the consideration of 
such a local under-density. As discussed in Section 12.41 
we note that the application of this method has the ef- 
fect of "correcting" the X— band galaxy counts from this 
sample (shown in Figure [6|), such that they agree with 
previously published results from the literature. 

While this method is rather crude, it allows for a means 
of estimating the light lost due to aperture clipping. We 
found that the median clipped aperture correction was 
~ 0.1 magnitudes. We applied this method to all galax- 
ies with clipped apertures down to 2 magnitudes below 
our selection limit for this study {K < 16.3). We found 
that this correction only increased the total number of 
galaxies in the sample by ~ 0.25% (via moving galax- 
ies from fainter to brighter than the magnitude selection 
limit after the aperture correction). 

2.3. Area Estimation 

Given that we wish to investigate the surface density 
and volume density of various populations of galaxies, we 
need to know the area on the sky (and hence volume) of 
our survey. Due to the complex shape of the gaps and 
holes in this survey, a "by eye" determination of area 
would be tedious and unreliable. Instead, we employ 
a "counts-in-cells" method of determining area coverage. 
To do this, we first find the distance to the nearest neigh- 
bor on the sky for all objects in our catalog. Given the 
typical (median) nearest neighbor distance of (5 « 150", 
we divide the entire area surveyed into square cells S x S 
in size. 

From the center of each cell we tested a range of radii 
to search for nearby catalog objects. If catalog objects 
were detected within the search radius, then the cell was 
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Figure 6. X— band galaxy counts as a function of apparent mag- 
nitude for th is s tudy compared with average Xj— band galaxy 
counts fromTCcenan ct al. (2010b). The counts have been divided 
through by a normalized Euclidean model of slope a = 0.6 to ex- 
pand the ordinate (assuming the bright galaxy counts take the 
form N(m) = A X 1 0""^, where A is a constant). Counts from 
IKeenan et al.l lj2010bl ') are shown in red asterisks. Raw counts from 
this study are shown as black diamonds. Counts from this study af- 
ter the aperture clipping correction was applied are shown in blue 
diamonds. These results indicate that the shape and normaliza- 
tion of X— band galaxy counts in this study are in good agreement 
with previous studies from the literature. This suggests that our 
area calculation is accurate, our star-galaxy separation is robust, 
and our ape rture corre ction is having the desired effect. Note: we 
shift the .Keenan et aP (2010b) Xs— band counts by -1-0.15 in mag- 
nitude to adjust for the typical magnitude difference for galaxies 
b etween Kg and the UK IDSS ii"— band filter following the results 
of lKeenan eFaLl l l2010bl ). 



counted as surveyed area. We then simply count up 
the number of surveyed cells to compute the total area 
surveyed. We determined an optimal nearest-neighbor 
search radius of 300". Using this method we compute a 
total area for this work of 536.9 deg^ on the sky. Varying 
our cell size and search radius by a factor of 2 was found 
to change our area estimate by only ~ 2%. Thus, we 
estimate that our error in total area calculation is ~ 2%. 
The exact estimate of the area is not critical to the re- 
sults presented below, as it only serves to set the overall 
normalization. 

2.4. Catalog Checks and Basic Analyses 

The catalog generation methods described thus far 
yielded a catalog of 35, 342 objects. Here we describe 
some of the basic consistency checks that we performed 
to verify that the final catalog was indeed what we in- 
tended it to be, namely, a K—hand selected sample 
{K < 16.3) of all galaxies with NIR and optical pho- 
tometry and spectroscopy where it existed. 

As a fundamental check that we indeed have generated 
a magnitude limited sample of galaxies that is neither se- 
riously contaminated by remaining stars nor plagued by 
the unintentional removal of a large number of galax- 
ies in the catalog cuts, we turn to the galaxy counts as 
a function of apparent magn itude in our fi e ld. We com- 
pare our counts with those of lKeenan et al.l (j2010br ). who 
recently combined their deep wide-field NIR counts with 
data from the literature to come up with the best current 
estimate of average galaxy counts over a wide range in 
apparent magnitude. 
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The results of this comparison are shown in Figure [HI 
Black data points show the counts in our sample before 
the aperture clipping correction was applied. The aper- 
ture corrected counts are shown in blue, and we note that 
the counts for the aperture corrected data agree well with 
previous counts from the literature. These results sug- 
gest that our area estimation is accurate, our star-galaxy 
separation is robust, and our aperture correction method 
is appropriate. 

2.5. Determination of Absolute Magnitudes 

Ultimately, we wish to make a comparison of the rest- 
frame NIR luminosity density as a function of redshift. 
To do this, we need to adjust the observed apparent mag- 
nitudes in our sample by a distance modulus (DM), a 
i^(z)— correction to correct for bandpass shifting, and an 
evolution correction, E{z), such that the absolute mag- 
nitudes used for constructing luminosity functions are 
given as: 



L 



^Q-0.i{M-M*) 



(3) 



M = m- DM{z) - K{z) + E{z). 



(1) 



At low redshifts in the NIR, J C(z)— corrections are n early 
independent of galaxy type (jMannucci et al.ll2001[ ). al- 
lowing for this magnitude correction to be made without 
considering the type distribution of the galaxy sample. 
Comb ining SDSS and UKIDSS data, iChilingarian et all 
(J201C0 showed that, at z < 0.5, accurate iiT— corrections 
can be calculated using low-order polynomials with in- 
put parameters of only the redshift and one observed 
NIR color. They have provided a iiT— correction calcula- 
tor packag43, which we used t o compute JiT— corr e ctions 
for galaxies in our sample. IChilingarian et all ()201Clf ) 
compared their iiT— correction calculator output values 
with those obtained via the more rigorous spectral energy 
distri bution fitting methods of Blanton fc RoweisI (|2007[ ) 
and Fioc & Rocc a-Volmerangel (J1997D and concluded 
that the magnitude errors associated with X— corrections 
derived using their algorithm should be < 0.1 magni- 
tudes. 

Evolution of the rest-frame NIR light from galaxies is 
expect ed to be significantly weaker than in optical band- 
passes (jBlanton et al.|[2"003[ ) , but it is an effect that must 
be accounted for when comparing galaxy luminosities at 
different redshifts. A commonly assumed form of the 
evolution corr ection is E(z) = Qz , where Q is a posi- 
tive constant. IBlanton et al.l (j2003f) showed that in the 
NIR, Q — 1 agrees well with stellar population synthesis 
models. Thus, for this study, we adopt Q = 1, such that 
E{z) — z. We further discuss this evolution correction, 
and its associated uncertainties, in Section 13.4.11 

3. THE E'-BAND GALAXY LUMINOSITY FUNCTION 

A number of different methods exist for estimating the 
galaxy LF. For a n excellent revi ew of the subject, we 
refer the reader to lJohnstoiil ()2011[ ). Th e most coinmonly 
assumed form of the LF is that of the iSchechten ()1976l ) 
function 

«I,.I..-(A)°„p(^)f, (2, 

which may be written in terms of absolute magnitudes 
using 

^ http://kcor.sai.msu.ru/ 



givmg 



$(Af ) = O.41n(lO)0* 



^qOA(M' -M) 



(Q + l) 



exp(100-4(M*-M)) 



(4) 



The Schechter function parameter L^ (or M*) repre- 
sents the luminosity of galaxies at the knee of the LF. 
(f>* determines the number density of L^ galaxies, and a 
is the faint-end slope. While the Schechter function has 
been demonstr ated to provide a less-than-perfect fit to 
real data (e.g., iJones et"al]|2006[) . it can provide a rea- 
sonably good fit and is the most widely used functional 
form for fitting the LF, making it the most useful form 
to consider when comparing with other studies from the 
literature. 

To fit Schechter functions to observed data, a va- 
riety of methods h ave been used in the past. In 
iKeenan et al.l (|2012f ). we compared four different LF 
estimators (l/Knax,C-, STY, and SWML) in the de- 
termination of NIR LFs. We found that the STY 
(iSandage et al.lll9790 and SW ML (Step- Wise Maximum 
Likelihood. lEfstathiou et aT]ll988.1 methods yielded sim- 
ilar resu lts in the determ ination of M* and a, while 
the C- (|Lvnden-Belll[l97l and 1/Knax (|Schmidtlll96"8h 
methods tended to underestini ate the faint end slope 
(see also IPage fc Carreral 120001 ) . Of these four meth- 
ods, l/14nax is the only one that p rovides the normaliza- 
tion {ip*). In IKeenan et al.l (j2012i ). we tested the 1/Vmax 
meth od alongside thr e e oth er normalization estimators 
from iDavis fc Huchral ()1982f ) . We found all four of these 
estimators yielded consistent results, a confirmation of 
the same result found by WiUmer (1997) using simulated 
data. 

Given these analyses, in IKeenan et al.l (j2012[) . we set- 
tled on a hybrid method to estimate the LF by first using 
STY to calculate Af * and a, and then using l/Knax (with 
M* and a fixed) to determine the normalization. Here 
we use this same hybrid method in the determination 
of the iC— band LF. To correct for spectroscopic incom- 
pleteness, we use a simple scheme in which each galaxy 
counted in the 1/T4nax procedure is weighted by a factor 
of 1/C{m), where C{m) is the fractional completeness as 
a function of apparent magnitude. 

3.1. The Assumption of a Constant LF Shape 

In what follows, we make the critical assumption that 
the z = shape parameters of the iiT— band LF (Af * 
and a) are not changing as a function of environment or 
distance from us. We require this assumption to facilitate 
the measurement of the JiT— band luminosity density as 
a function of redshift over this range. This is due to the 
fact that, in this sample, we have little or no information 
about the faint end slope of the LF at higher redshifts 
due to the magnitude limit of the survey, and we have 
corrupted information about the bright end of the LF 
at lower redshifts due to the Petrosian aperture clipping 
issue described in Section [^21 

We believe the assumption of a constant LF shape 
is reasonable, given that the iir(z)— corrections are es- 
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Figure 7. The UKIDSS-LAS A'— band luminosity function over 
the range 0.005 < z < 0.2. The red circles show the 1/Vmax LF 
estimate for this redshift range. Error bars show Poisson counting 
errors. The black curve shows the best fit normalization, 0* = 
0.015, to l/Vmax data after having determined the values of M* = 
-21.37 and a = -0.89 using the STY method. The value of x^^j = 
4.64 demonstrates the relatively poor fit of the Schechter function 
to the data. 



sentially independent of galaxy type and the E{z) cor- 
rectio ns are quite modest (see Sect ion 12.51) . Fu rther- 
more.JDe Propris fc ChristleinI (|2009( ) and Capozz i et al.l 
(|2012[ ) have found that the shape of the NIR LF is not 
significantly different for field and cluster galaxies. Thus, 
in the following sections, we consider the K—hand LF for 
all galaxies in our sample under the assumption that the 
shape of the z = LF is not changing significantly as a 
function of distance from us or environment. 

The assumption of a constant LF shape allows us to 
derive Af * using the higher redshift galaxies in the sample 
and then to determine the faint end slope, a, using the 
lower redshift galaxies. Thus, we used the STY method 
on the sample at redshifts 0.09 < z < 0.2 to determine 
M* with a fixed. We then used STY on the sample 
at 0.005 < z < 0.07 to determine a with M* fixed at 
the value determined using the higher redshift sample. 
We iterated this process until we reached convergence 
in the values of Af * and a. With M* and a fixed, we 
then investigated the changing normalization (and hence 
changing luminosity density) as a function of redshift. 

3.2. The K—band luminosity density at 0.005 < z < 0.2 

In Figure [3 we show the K—ha.nd LF for all galaxies 
in the redshift range 0.005 < z < 0.2. The red circles 
show the 1/Vrnax estimate of the LF, and the black curve 
shows a fit to these data having determined Af * and a 
using STY and then fitting 0* to the l/Vmax data. We 
note the poor fit to the data (Xred = ^■^)- ^ poor fit of 
the Schechter function to LF estimate s is not uncommon 
in studies from the liter ature (see e.g.. lLavaux fc HudsonI 
120111 : iSmith et al.ll2009[ ) , and authors normally attribute 
this to the Schechter function being a poor model of the 
LF. However, in the analysis that follows we demonstrate 
that much of this fitting error may arise from density 
variations along the line of sight, and that in smaller red- 
shift bins, the Schechter function provides a reasonable 
fit to the data. 

It is worth noting at this point that a considerable por- 
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Figure 8. The UKIDSS-LAS ii'-band luminosity function split 
into three redshift ranges: 0.005 < z < 0.07 (red), 0.07 < z < 0.09 
(green), and 0.09 < ^ < 0.15 (blue). We separate the redshift range 
0.07 < z < 0.09 to demonstrate that the excess at z > 0.09 is not 
due to the Sloan Great Wall or the other over-densities we observe 
at higher declination in this redshift range. Here we have deter- 
mined the LF shape parameters (a = —1.02 and Af* = —21.38) 
itcratively using STY on b oth the high and low-redshift samples 
as described in Section |3. II We then fit <j>* in each redshift range 
with a and A/* fixed. We list the X^cd values for each redshift 
range in the figure, and note that relatively good fits can be ob- 
tained using fixed LF shape parameters and letting the normal- 
ization (</<*) vary as a function of redshift. We attribute the poor 
fit at 0.07 < z < 0.09 to the large inhomogeneities in this range. 
However, this redshift bin was not use in the determination of the 
LF shape. As a result of these analyses, we find a value for </<* that 
is ~ 1.4 times higher at 0.09 < z < 0.15 than at 0.005 < z < 0.07, 
suggesting a relative luminosity density that is high by the same 
factor. Error bars represent Poisson counting errors. 



tion (^ 25%) of the solid angle subtended on the sky by 
this s ample is filled by the Sloan Great Wall (Go tt et al.l 
120051 ) at redshifts of 0.07 < z < 0.09. In our sample, 
^ 40% of the galaxies in this redshift range are part of 
this structure. The relative excess in the redshift dis- 
tribution due to the Sloan Great Wall can be seen as a 
peak at these redshifts in Figure [3l both in the overall 
redshift distribution and, more prominently, in subregion 
2, which is centered on the this structure. Such struc- 
tures will certainly be expected to cause some deformity 
in the LF, given that different luminosity ranges are pref- 
erentially sampled at different redshifts in an apparent 
magnitude limited survey. 

In the same sense, it is worth noting that in Figure [7l 
the faint end of the LF is predominantly coming from 
the low-redshift end of the sample, and the bright end of 
the LF is mainly coming from the higher redshifts. With 
this in mind, in Figure [51 we show the LF for low-redshift 
galaxies in the range 0.005 < z < 0.07 as red circles, the 
LF for the redshift range 0.07 < z < 0.09 as green circles, 
and the LF for galaxies at 0.09 < z < 0.2 as blue circles. 
The LF shape parameters, M* and a, are the same for 
all three LFs shown in this figure. We determined M* 
and a iteratively using the STY method by fitting M* on 
the high-redshift (0.09 < z < 0.2) sample with a fixed, 
then fitting a on the low-redshift (0.005 < z < 0.07) 
sample with M* fixed. We repeated this procedure until 
we reached convergence at values of Af * = —21.38 ± 0.04 
anda = -1.02 ±0.03. 

With M* and a fixed, we fit the normalization, 0*, 
to the LF estimate given by the l/14iax method. We 
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Figure 9. UKIDSS-LAS X— band LFs as a function of redsiiift. Here wehave fixed tfie LF shape parameters (a = —1.02 and M* = —21.38) 
to those derived from the iterative fitting method described in Section |3. 1 1 using the data presented in Figure |8] We then divided the sample 
into the six independent redshift bins shown here. In each case, the error bars show Poisson counting errors. We find that we can get 
a relatively good fit over all redshift bins (with the notable exception of the 0.07 < z < 0.09 bin containing the Sloan Great Wall) by 
simply letting the Schechter function normalization vary as a function of redshift. The normalization {(j>*) and number of galaxies used to 
determine the LF in each bin are listed in the plots, as well as the X?cd values for each case. 
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Figure 10. UKIDSS-LAS ii"— band luminosity density as a function of redshift. (a) Our measured ii"— band luminosity density for the full 
sample (red circles) versus different directions on the sky (green, blue, and orange circles corresponding to subregions 1, 2, and 3, respectively, 
in Figure [!}• (b) Our measured X— band luminosity density for the full sample (red circles) as a function of redshift compared with other 
studies from the literature. The luminosity density for studies from the literature have been recalculated here using a consistent Mq k 
and Schechter function integration limits in all cases. The error bars on the data from this study in panel (b) represent the 1 a statistical 
errors in luminosity density in each bin combined with the systematic errors due to cosmic variance estimated from panel (a) of this plot. 
The density contrast, (p)/po, is displayed on the right-hand vertical axis. The scale of the right-hand axis was established by performing 
an error-weighted least-squares fit (for the normalization only, not shape) of the void radial density profile of Boleiko & Sussman (20lJ) 
(gray solid curve) to all the luminosity density data in pan el (b). Th e dashed curve shows the void radial density profile of Alexande r et all 
l|2009). Both Alexander ct al. (2009) and Bolejko & Suss manI 1)20111 ) claim these density profiles can provide for good fits to the SNI a data 
without dark energy. The dash-dot curve shows the scale and amplitude of the "Hubble bubble" type perturbation that Marra et al] l|2013l ) 
would require to explain the discrepancy between local measurements of the Hubble constant and those inferred by Planck. While all of 
the models shown for comparison are essentially toy models, they demonstrate the rough size scale of the under-densitios considered in 
current minimal-void cosmological models. 
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Table 1 

NIR luminosity densities and Schechter functions parameters for this study and selections from the literature 





{z)- 


Redshift Range 


^gals 


M*^ 


a 


(/)*<^xl03 


iK,calc'* 


jK.pub 


This Study 


0.028 


0.005 <z< 0.04 


2,351 


-21.38 ±0.04 


-1.02 ±0.03 


10.9 ±0.5 


5.15±0.44(±0.31) 


— 


This Study 


0.055 


0.04 < 2 < 0.70 


5,185 


-21.38 ±0.04 


-1.02 ±0.03 


12.1 ±0.6 


5.73±0.41(±0.81) 


— 


This Study 


0.081 


0.07 < 2 < 0.09 


6,303 


-21.38 ±0.04 


-1.02 ±0.03 


20.1 ±0.8 


9.47±0.65(±0.97) 


— 


This Study 


0.099 


0.09 < 2 <0.11 


4,693 


-21.38 ±0.04 


-1.02 ±0.03 


16.0 ±0.4 


7.53±0.55(±0.75) 


— 


This Study 


0.119 


0.11 < 2 < 0.13 


4,107 


-21.38 ±0.04 


-1.02 ±0.03 


16.6 ±0.5 


7.84±0.58(±0.41) 


— 


This Study 


0.154 


0.13 < 2 < 0.20 


8,095 


-21.38 ±0.04 


-1.02 ±0.03 


16.5 ±0.4 


7.83±0.51(±0.29) 


— 


Keenan+ (2012) 


0.200 


0.10 < 2 < 0.30 


812 


-21.56 ±0.06 


-0.91 ±0.07 


14.6 ±1.7 


7.48 ±1.03 


— 


Lavaux+ (2011) 


0.030 


0.017 < 2 < 0.067 


60, 000 


-21.33 ±0.02 


-0.73 ±0.02 


11.1 ±0.2 


3.81 ±0.25 


(3.94 ± 0.02) 


HillH- (2010) 


0.076 


0.003 < 2 < 0.1 


1,785 


-21.66 ±0.10 


-0.96 ±0.06 


15.6 ±1.6 


7.65 ±2.16 


(6.98 ±1.62) 


Smith+ (2009) 


0.100 


0.01 < 2 < 0.3 


40,111 


-21.48 ±0.04 


-0.81 ±0.04 


16.6 ±0.8 


6.55 ±0.84 


(6.18 ±0.07) 


Jones+ (2006) 


0.054 


0.0025 < 2 < 0.15 


60, 869 


-21.98 ±0.03 


-1.16 ±0.04 


7.5 ±0.5 


5.80 ±0.85 


(5.80 ±0.50) 


Eke+ (2005) 


0.070 


0.005 < 2 < 0.12 


15,644 


-21.58 ±0.04 


-0.81 ±0.07 


14.3 ±0.8 


6.30 ±0.86 


(7.04 ± 0.23) 


Bell+ (2003) 


0.078 


0.0033 < 2 < 0.2 


6,282 


-21.44 ±0.05 


-0.77 ±0.04 


14.3 ±0.7 


5.48 ±0.73 


(5.80 ± 1.80) 


Huang+ (2003) 


0.138 


0.005 < 2 < 0.35 


1,056 


-21.85 ±0.08 


-1.37 ±0.10 


13.0 ±3.0 


11.32 ±2.56 


— 


Feulner+ (2003) 


0.200 


0.10 < 2 < 0.30 


210 


-21.94 ±0.24 


-1.10 ±0.10 


11.1 ±1.2 


7.91 ±2.00 


— 


Cole+ (2001) 


0.060 


0.023 < 2 < 0.12 


5,683 


-21.59 ±0.03 


-0.96 ±0.06 


10.8 ±1.6 


5.10 ±1.29 


(5.74 ±0.86) 


Kochanek+ (2001) 


0.023 


0.007 < 2 < 0.1 


3,878 


-21.54 ±0.05 


-1.09 ±0.06 


11.6 ±1.0 


6.87 ±1.62 


(7.14 ±0.75) 
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find that the LFs for the high-redshift and low-redshift 
samples shown in Figure [8] can be reasonably well fit by 
the same M* and a parameters, while only allowing the 
normalization to change (we attribute the comparably 
poor fit at 0.07 < z < 0.09 to the large inhomogeneities 
in this redshift bin). This indicates that the assumption 
of a constant LF shape as a function of redshift and en- 
vironment appears to be valid. As shown in Figure [8l 
we find a normalization at 0.09 < z < 0.2 that is ^ 1.4 
times higher than that at 0.005 < z < 0.07. Thus, if 
the LF shape is indeed constant for this sample, then 
the luminosity density at z > 0.09 appears to be ~ 1.4 
times higher than that at z < 0.07, even when known 
over-densities at 0.07 < z < 0.09 are excluded. 

Next, we investigated the ii'— band luminosity density 
as a function of redshift in a series of narrower redshift 
bins over the range 0.005 < z < 0.2, as labeled in each 
panel of Figure [9] We fixed the LF shape parameters to 
the values derived above (a = -1.02 and M* = -21.38) 
and fit the normalization to the 1/Vmax results in each 
redshift bin. Red circles show the 1/Vinax results in each 
case, and the black solid curves represent the Schechter 
function fits. Again, we find that the data in all of the 
redshift bins are reasonably well fit with the same M* 
and a parameters (x^ed "^ ^ except at 0.07 < z < 0.09), 
suggesting that the assumption of a constant LF shape 
is valid. We also note in each panel the measured value 
for (/)* and the number of galaxies used to generate the 
LF. 

3.3. Comparison With Studies From the Literature 

In Figure [ini we use the results from Figure IH] to calcu- 
late the if— band luminosity density as a function of red- 
shift (by integrating the fitted Schechter functions). In 
Figure llOb , we show how these results vary as a function 



of position on the sky by comparing the average result for 
the full sample (red circles) with the results from each of 
the three subsamples (green, blue, and orange circles) of 
galaxies from the subregions shown in Figure [1] Each of 
these subsamples contains roughly a third of the original 
sample of galaxies. While we find stark differences be- 
tween subsamples in the measured luminosity density in 
some redshift bins, the overall trend toward a rising lu- 
minosity density with increasing redshift appears present 
in all cases. We use the rms variation in luminosity den- 
sity in each redshift bin as an estimate of the systematic 
error due to cosmic variance in our measurement. These 
systematics are refiected in the larger error bars for the 
total sample in Figure fTOb. 

In Figure [TOb . we compare our results with other stud- 
ies from the literature, where red circles again show our 
results for the entire sample, and studies from the litera- 
ture correspond to the symbols listed in the plot legend. 
The luminosity densities presented in this figure for stud- 
ies from the literature have been recalculated here using 
the same value for Mq^k and Schechter function inte- 
gration limits in all cases for consistency. All the data 
displayed here are listed in Table [T] The majority of 
studies at z < 0.1 have used photometr y from th e Two 
Micron Ah Sky Survey (2MASS, Jarr ett et al.l [20001. 
These studies have typically used the 2MASS ifs— band 
(2.12/im) Kron or "total" magnitudes which are gener- 
ally found to be ^ 0.15 mags brighter than UKIDSS 
if —band Petrosian magn itudes (see e.g., iKeenan et alj 
I2010at iSmithet al.ll2009[ ). Thus, in making this com- 
parison, we have adjusted the values for M* to be 0.15 
mags brighter for o ur study and other studies d erived 
from UKIDSS data (|Hin et al.ll20"Tfll : ISmith etHI [20091. 
Making this adjustment does not change the results pre- 
sented here, but is done for the sake of consistency. In 
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general, we find excellent agreement with all previously 
published results in the -ftT— band, and we confirm the 
tentative result presented in iKeenan et al.l ()2012[ ) of a 
rising luminosity density from z ~ 0.05 to z ~ 0.1, which 
appears to remain higher than that measured locally out 
to z~ 0.2. 

The relative density contrast, (p)//5o, is displayed on 
the right-hand vertical axis in Figure [TOl The scale 
of this axis was established by performing an error- 
weighted least-squares fit of the void radial density pro- 
file oflBolciko & Sussman (2011) (solid curve, fitting only 
the normalization only, not the shape) to all the luminos- 
ity density data in panel (b). As mentioned in Section [U 
this means we are assuming a linear bias parameter of 
6 = 1, or more specifically that ii'— band luminosity den- 
sity is an unbiased tracer of the underlying dark matter 
density. We believe this is a reasonable assumption given 
previous results from observations and simulations (e.g., 
ikaller et al1[200llAM^lo et al.llMl . 

The dashed cur ve in Figure [TOl s hows the void radial 
density profile of [Alexander et al.l ()2009[) , also normal- 
ized to a density contrast of unity at high redshift. Both 
[Alexander et alJ (|2009) and Bolcjko & Sussman (20lil) 
claim these density profiles allow for a good fit to the 
SNIa data without dark energy. The dash-dot curve 
shows the scale and a mplitude of the "H ubble bubble" 
type perturbation that iMarra et al.l ()2013[ ) would require 
to explain the discrepancy between local measurements 
of the Hubble constant and those inferred by Planck. 
While all three of these models are essentially toy mod- 
els, they demonstrate the rough size scale of the under- 
densities considered in current minimal- void cosmologi- 
cal models. 

In general, the data from this study and those from the 
literature seem to indicate a rising luminosity density out 
to z ~ 0.1. The not able exception to this t rend would be 
the measurement of lKochanek et ahl (|200lD . However, we 
note that the relat i vely lo w luminosity density found by 
iLavaux fc Hudson! ()20lH ) is also derived from 2MASS, 
but with much improved photometry, and includes a 
sample more than an order of magnitu de larger, which 
itself contains the original sample from iKochanek et al.l 
((2001 . , , , 

The stu dies of ICole et~all (12001. IBeU et al.l (l200l . 
lEke et al.l (|2005[ ). and iJones et al.l (|2006[ ) are also based 
on 2MASS photometry, and all agree well on a signifi- 
cantly hi gher value for the lumino sity density than that 
found bv lLavaux fc HudsonI (|2011| ). However, all four of 
these studies considered deeper redshift samples extend- 
ing to z > 0.1. Qualitatively speaking, in all four cases 
a tension similar to that which we show in Figure [71 be- 
tween the observed LF and the best Schechter function 
fit, appears present. Thi s phenomenon is se en most dra- 
matically in the work of IJones et al.l (|2006l ). where they 
display the residuals of their Schechter fits to the LF 
showing that the faint end is overestimated by the fit, 
while there appears to be an excess of galaxies around 
M ~ M*. We would argue, based on our own study of 
the LF, that much of the irregularity in the LF shape 
found by these studies may be due to inhomogeneity 
along the line of sight, and more specifically, to a higher 
space density of galaxies at the high-redshift end of the 
sam ples. 

.Hill et all (|2010D use UKIDSS data to measure the 



ii'— band luminosity density in the field targeted by the 
Millennium Galaxy Catalog Survey (Dec = 0, lO'' < 
RA < IS^. lLiskeet al.ll2003| ). As noted in Section [O 
this r egion contains the Sloan Great Wall. Hill et aQ 
()2010t ) correct for this over-density by adjusting their 
normalization down by ~ 20%. Undoing this correction 
would bring their measurement into good agreement with 
our result at roug hly the same redshift. 

The sample of iSmith et al.l ([2009) is the most simi- 
lar to our own. They select ^ 40, 000 galaxies from the 
UKIDSS-LAS with redshifts from the SDSS. They mea- 
sure the LF for the entire sample from 0.01 < z < 0.3. 
The find similar irregularities in the LF shape to those 
we present in Figure [T) They note that when they divide 
their sample into three redshift bins they see a rising 
LF normalization with increasing redshift, and point out 
that even doubling their evolution correction does not 
resolve this issue. 

The samp le sizes of iHuang et ahl (|2003D . iFeulner et al.l 
()2003() , and IKeenan et al.l ()2012D are generally too small 
to be considered robust to cosmic varianc e. However, it 
is worth notin g that the measu rements of IFeulner et al.l 
(|2b03) and Keenan et al.l (|2012[ ) may be considered con- 
servative underestimates of the true luminosity density 
because these studies avoided known over-densities, such 
as galaxy clusters, in the redshift ranges sampled. 

We conclude that if the observed trend in i^T— band lu- 
minosity density as a function of redshift is indicative of 
a similar trend in the underlying total mass density, then 
the local universe may be under-dense on a scale and am- 
plitude sufficient to introduce significant biases into lo- 
cal measurements of cosmological observables. Leaving 
aside considerations of whether or not such an unusual 
local structure coul d obviate the need for dark energy, 
given the analysis of IMarra et al.l (|2013l ). it appears that 
the observed under-density is roughly the right scale and 
amplitude to explain the apparent tension between local 
measurements of the Hubble constant {Hq = 73.8 ± 2.4 
km s~^ Mpc~^, iRiess et all 120 lit ) and the recent re- 
sults from Planck (Hn = 67.3 ± 1.2 km s'^ Mpc-\ 
iPlanck Collaboration et al.ll2013[ ). 

3.4. Possible Sources of Bias and Error in This 
Measurement 

We have done everything possible in this study to make 
an unbiased measurement of the /C— band luminosity 
density as a function of redshift in the nearby universe. 
However, we are pushing the data by making several as- 
sumptions that allow us to probe a wider redshift range 
than that for which we have a statistically robust mea- 
surement of the LF over the entire range of absolute mag- 
nitudes. Here we discuss the possible biases and errors 
in our measurement associated with these assumptions. 

3.4.1. K(z) and E(z) corrections 

The most important assumption we make in this study 
is that the 2 = shape of the K—hand LF is not chang- 
ing significantly as a function of distance from us. This, 
in turn, relies on the assumption that we are making 
the appropriate K{z) and E{z) corrections to calculate 
the z = Q absolute magnitudes of galaxies at any given 
redshift. To investigate these assumptions in detail, we 
explored variations in the K{z) and E{z) corrections to 
see how they affect our result. 
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If we simply omit the i^(z)— corrections in our cal- 
culations, we find that the excess luminosity density 
at z > 0.1 shown in Figure [10] increases significantly 
(due to the fact that the ii'— correction is negative). 
We compared the /^—co rrection calculator outputs of 
Chilingarian et aTl (120101) to the empirical estimates by 
Mannucci et al.l (|2001[ ) and found qualit ative agreement, 
but o n average the iiT— corrections of IMannucci et al.l 
()2001| ) are ^ 0.1 mags more negative than the calcu- 
lator outputs. We ran all of our analyses using the aver- 
age A'— corrections of Mannucci et al. (2001) and found 
that the excess luminosity density at ^ > 0.1 is reduced 
slightly, but not significantly. In general, we expect the 
calculator outputs to be more robust estimates of the 
true if— corrections as the y were derived u s ing a much 
larger sample than that of IMannucci et aLl ([2001'), and 
the sample was drawn from the UKIDSS and SDSS sur- 
veys, which we use for this study. 

The E{z) corrections we use are of the form E{z) = 
Qz, as described in Section [231 and also act to slightly 
reduce the observed excess in luminosity density. To ex- 
plore the possibility of underestimated evolution, we in- 
creased the Q parameter arbitrarily to investigate the 
effects of stronger evolution. We found we needed to in- 
crease the evolution correction to Q = 4 to make the 
average luminosity density at z > 0.1 roughly equal to 
that at z = 0.05. This would imply that galaxies are 
dimming by ^ 30% in the rest-frame K—hand over the 
last 1 — 2 Gyrs. This amount of evolution would require 
galaxies to form at 2; <C 1, which is ruled out by galaxy 
star formation histories. Thus, we do not expect under- 
estimates of the E{z) corrections to be a possible source 
of the excess luminosity density at z > 0.1. 

3.4.2. The Faint-end Slope of the Luminosity Function 

At z > 0.1 we are not making a robust estimate of the 
faint end slope (a) of the LF. We examined the possi- 
bility that a less negative slope than that measured at 
z < 0.07 (i.e., a < —1.02) at z > 0.1 could account 
for the observed excess in luminosity density. This is, 
essentially, another means of exploring the possible ef- 
fects of unanticipated evolution in the LF as a function 
of redshift. We found that a slope of a w —0.65 im- 
posed at z > 0.1 (while leaving a = —1.02 at z < 0.07) 
was sufficient to suppress the observed excess at z > 0.1. 
However, this would not explain the rising normalization 
with increasing redshift. In addition, the vast majority of 
studies of the NIR LF from the literat ure measure a to 
be ~ 0.9 - 1.1 for < z < 0.3 Ce.g.. ICole et al.l[200l 

et al.l 
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Thus, a 



-0.65 at z > 0.1 



would not only represent some kind of extreme evolu- 
tion in the LF but also be at odds with results from the 
literature. Therefore, we do not consider this to be a 
reasonable candidate for the observed excess luminosity 
density at z > 0.1. 

3.4.3. Spectroscopic Selection Bias 

The spectroscopy for this study comes mainly from the 
SDSS and is supplemented by other published redshift 
surveys. Targets for SDSS spectroscopy were chosen us- 
ing an _R— band selection, and, in general, the targets 
for other redshift surveys from the literature were also 
selected optically. Thus, there could exist a bias in the 



spectroscopic redshift distribution for a AT- band selected 
sample for which spectroscopic targets have been selected 
in optical bands. At z ~ 0.05, the typical R— K color of 
galaxies is ~ 2, while at z ~ 0.15, i?— K ^ 1.25. Thus, 
there will be a bias against faint ii'— selected sources at 
low redshifts being chosen for spectroscopic follow-up in 
an _R— band target selection (relative to the same selec- 
tion at higher redshifts). In principle, such a selection 
bias could reduce the overall A'— band LF normalization 
at low redshifts. 

For this reason, we restrict this study to only areas 
of very high spectroscopic completeness. The minimum 
acceptable completeness for regions covered in this study 
is 90%, and, on average, the completeness is ~ 95% at 
ATab < 16.3. To estimate the limiting case of the R — 
K spectroscopic selection bias, we consider the extreme 
scenario in which all galaxies lacking spectroscopy are 
assumed to reside at z < 0.07. In this case, we estimate 
that the LF normalization at z < 0.07 could be increased 
by as much as ~ 20%. 

However, the median apparent AT— band magnitude of 
a galaxy at z < 0.07 is roughly 0.75 mags brighter than 
that of a galaxy at z > 0.09 in this study, so the spectro- 
scopic selection effect should be essentially cancelled out 
by the fact that nearby galaxies are intrinsically brighter 
on the sky, and we would expect any biases due to such 
selection effects at low redshifts to be much smaller than 
the 20% quoted above. 

If we assume that spectroscopic selection bias is not a 
problem in this study, then we can expand the area on 
the sky to include the entire region of overlap between 
the UKDISS-LAS and the SDSS. In doing so, the aver- 
age overall completeness of the sample drops to ^ 85%, 
but the volume is increased by roughly a factor of four 
(increase in area from ^ 500 deg^ to ^ 2000 deg^). We 
performed all of the analyses described above on this full 
sample of ^ 140, 000 galaxies and obtained very simi- 
lar results, namely, that the AT— band luminosity density 
at z > 0.1 appears to be ~ 1.5 times higher than that 
measured locally. 

3.4.4. Photometric Errors 



As noted in Section [521 galaxies which subtend a large 
solid angle on the sky may have their Petrosian fluxes un- 
derestimated due to the upper limit of 6" (corresponding 
to a 12" circular aperture radius) on the size of the Pet- 
rosian aperture in the WSA pipeline. As demonstrated 
in Figure [11 this issue presents a problem at low redshifts 
(z < 0.1) where a significant fraction of objects have had 
their Petrosian apertures clipped, and thus their fluxes 
underestimated. 

In Section 12.21 we describe a method to estimate a 
flux correction for galaxies where the Petrosian aperture 
was clipped by using a range of circular aperture magni- 
tudes to measure a light curve for each galaxy, and then 
extrapolating to the z— band Petrosian aperture radius 
given in the SDSS. We show in Figure [6] that making 
this correction brings the galaxy counts for our sample 
into agreement with the expected distribution taken from 
all-sky galaxy counts in the AT- band. We explored the 
effect of arbitrarily increasing this correction and found 
that a larger correction (by a factor of ~ 2) had the ef- 
fect of overpopulating the bright end of the galaxy counts 
distribution. 
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We found that applying the aperture clipping correc- 
tion factor before making the initial apparent magnitude 
selection in the X— band only increased the sample size 
by ~ 0.25%. Before making this aperture correction, 
we found that the bright end of the LF at low redshifts 
(z < 0.05) appeared to fall off more steeply than expected 
(given the distribution seen at slightly higher redshifts). 
Upon making the correction, the bright end of the LF 
at low redshifts was moved to brighter magnitudes and 
toward better agreement with higher redshift results. 

We note, however, that in the lowest redshift bin pre- 
sented in Figure [HI the bright end of the LF still appears 
to fall off more steeply than at higher redshifts. We found 
that arbitrarily increasing the aperture correction in this 
redshift bin appears to bring the bright end of the LF 
into better agreement with the LFs in higher redshift 
bins. Thus, we conclude that we are probably under- 
estimating the aperture clipping correction at the very 
lowest redshifts. However, the aperture correction itself 
(even when increased arbitrarily) does almost nothing to 
increase the normalization of the LF. This is because in- 
trinsically bright (large) galaxies are those most affected 
by aperture clipping, and correcting for this effect only 
shifts the bright end of the LF along the abscissa. 

Thus, while this method represents a rather crude lost 
light correction for any individual galaxy, we believe that 
it provides a statistically sound means of recovering the 
galaxy luminosity distribution. Furthermore, we find 
that the relative distribution of K—hand luminosity den- 
sity as a function of redshift presented in Figure ITOl is not 
significantly altered by making this correction (or dou- 
bling it for that matter). 

3.4.5. Cosmic Variance 

Systematic bias due to cosmic variance is a possible 
source of error in any measurement made over a volume 
that is insufficient to average over large-scale structure 
in the universe. As noted in Section [H the largest struc- 
tures observed in the local universe appear to be of order 
the size of the largest volumes surveyed to date. Thus, 
it remains unclear what the upper limit on the size of 
structure is and what volume constitutes a representa- 
tive sample of the universe. 

To address the issue of cosmic variance, we divided 
our sample into three subsamples denoted by the differ- 
ent regions on the sky shown in Figure [TJ As shown in 
Figure 1101 our measured luminosity density in some red- 
shift bins varied dramatically among the three subsam- 
ples, but, in general, the overall trend toward a rising 
luminosity density with increasing redshift is present in 
all three regions on the sky. 

In the highest redshift bins the results from the differ- 
ent subsamples appear to be converging, suggesting that 
some kind of 'scale of homogeneity' has been reached, but 
the substantial increase in luminosity density compared 
to the lowest redshift bins implies that structure may ex- 
ist on scales larger than the survey volume. We estimate 
the systematics due to cosmic variance in each of our 
six redshift bins given the rms scatter between measure- 
ments from the three subregions. We present this error 
estimate in Table [T] and graphically in Figure [TUb. 

As a final check that our methods are sound, and to 
investigate further the expected effects of cosmic vari- 
ance, we performed all the same analyses described above 



to estimate the relative variation in ii'— band luminos- 
ity density as a function of redshift using mock catalogs 
from the Millennium database onlinj^. The Millennium 
database contai ns a variety of mock catalogs based on 
the Millennium ('Soringel et al. 2005) and Millennium-II 
(Bovlan-Kolchi n et al. 2009) simulations. We used the 
"Blaizot2006" all-sky catalogs, which wer e generated us- 
ing the MoMaF code ()Blaizot et al.ll2005[) . as well as the 
"Henriques2012" all-sky ca talogs created us i ng th e vi- 
sual observing algorithm of Hcnrioucs ct al. ( 20121 ) and 
the semi-analytical model of Guo ot al. (2011). 

From both catalogs, we retrieved all-sky iiT— selected 
catalogs of galaxies at if < 16.3. We divided each all- 
sky catalog into 16 different regions to be considered 
independently. While each of these regions of ^ 2500 
deg^ is much larger than our observed sample, we ex- 
pect the effects of cosmic variance to be similar, be- 
cause our observed sample spans many widely separated 
fields on the sky. While it is difficult to quantify exactly 
how cosmic variance shoul d scale in such a comparison, 
[Driver fc RobothamI ()2010fl show that for a sample made 
up of N independent fields, cosmic variance should be re- 
duced by a factor of l/vTV. 

We then computed the iiT— band luminosity density as 
a function of redshift for each of the 16 regions in both 
mock catalogs using the same methods as were used with 
the observed data, namely, applying the same K{z) and 
E(z) corrections and estimating the LF in the same way. 
In Figure [m we show the results of this exercise, where 
blue circles represent the Henriques2012 mock catalog, 
and red circles represent the Blaizot2006 catalogs. The 
data points show the average result in each redshift bin, 
and all results are normalized to the average in the high- 
est redshift bin. The error bars show the standard devia- 
tion of the distribution of measured luminosity densities 
in each redshift bin for the 16 regions considered in each 
catalog. 

Thus, we find that a typical deviation of ^ 10 — 20% 
could be expected in this measurement between the low- 
est and highest redshift bins. While this result could 
be further quantified by including additional or indepen- 
dent realizations of mock catalogs, our interest here is 
simply to show that when we apply our methods to sim- 
ulated catalogs, we recover the expected result that the 
K—hand luminosity density is not changing dramatically 
as a function of redshift. 

4. SUMMARY 

We have presented a study of the if —band luminosity 
density in the redshift range 0.005 < z < 0.2. This study 
is based on a if— band selection of galaxies taken from 
the UKIDSS-LAS, where spectroscopy from the SDSS 
and other redshift surveys provides for > 500 deg^ on 
the sky that is > 90% spectroscopically complete to 
Kab < 16.3. The primary motivation for this study is to 
explore the low-redshift stellar mass density distribution 
of the universe to test for the possibility that the nearby 
universe is under-dense on a scale and amplitude suffi- 
cient to introduce significant biases into locally measured 
cosmological observables. This work represents the first 
time that the NIR luminosity density has been measured 

^'^ http;//gavo. mpa-garching.mpg.de/MyMillennium/ 
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Figure 11. Normalized luminosity density in the .ff— band de- 
rived from two different mock catalogs using all the same methods 
as were used on the observed data. In each case, we divided the 
all sky catalog for the K < 16.3 sample into 16 sight lines, for 
which we calculated the luminosity density in each of the six red- 
shift bins individually. The data points show the average result 
in each redshift bin for the Blaizot2006 catalogs (red circles) and 
the Henriques2012 catalogs (blue circles). The error bars show the 
standard deviation of the results over the 16 sight lines in each 
redshift bin for both catalogs. The results are shown normalized 
to the average in the highest redshift bin. 



the ne arby universe. Given the analysis of iMarra et aLl 
(|2013D . we can say that such a scenario would be suf- 
ficient to explain the apparent tension between direct 
measurements of the Hubble constant and those recently 
published by the Planck collaboration. 

As NIR luminosity is a good tracer of stellar mass, 
and stellar mass density should trace the underlying dark 
matter distribution, we c ompare our results w i th the void 
radial density profiles o f [Boleiko fc SussmanI ()201lD and 
[Alexander et al.l (|2009t ). which they claim could act to 
produce the apparent acceleration of the expansion of the 
universe observed via type la supernovae. While these 
void models are simplistic, they give a rough idea of the 
minimal scale and amplitude required for a local under- 
density to introduce significant biases in locally measured 
cosmological observables, such as the expansion rate. We 
conclude that local measurements of the NIR luminosity 
density could be described as consistent with cosmolog- 
ical models that invoke a large local under-density on 
radial scales of a few hundred Mpc to explain the appar- 
ent acceleration of the expansion of the universe. 



in the nearby universe with consistent methods and pho- 
tometry over a wide range in redshift. 

We find that we can achieve good fits to the iiT— band 
galaxy LF both at low redshifts and at higher redshifts 
using the same LF shape parameters, M* and a. We find 
values for these LF shape parameters of M* = — 21.38 ± 
0.04 and a = — 1.02± 0.03 by iteratively fitting a on the 
low-redshift data while holding M* fixed, then fitting M* 
on the higher-redshift data while holding a fixed until 
convergence is achieved. 

Next, we consider six redshift bins over the range 
0.005 < z < 0.2 to explore the X— band luminosity den- 
sity as a function of redshift in detail. To do this, we fix 
M* and a and only fit the normalization, 0* , in each bin. 
We find excellent agreement with all previous measure- 
ments of the iC— band LF from the litera ture. We con- 
firm t he tentative detection, presented in iKeenan et al.l 
(|2012l) , of a rising NIR luminosity density over the range 
0.005 < z < 0.1 and a luminosity density at z ^ 0.2 that 
is ^ 1.5 times higher than that measured locally. 

We also divide our sample into three subsamples to 
consider the issue of cosmic variance. We find that, while 
individual measurements of the luminosity density in a 
given redshift bin may differ by a factor of two between 
subregions, the overall trend with redshift is consistent 
between the subregions. At the highest redshifts consid- 
ered in this sample, the measurement in the three sub- 
samples appears to be converging on the result that the 
luminosity density ai z > 0.1 appears to be '-^ 1.5 times 
higher than that measured locally. In particular, it is 
worth noting that in subregions 1 and 2, which represent 
opposite directions on the sky, the measurements agree 
very well on the luminosity density in both the lowest and 
highest redshift bins. This relative agreement between 
widely separated sight lines, combined with the fact that 
the iiT— band galaxy counts from this study agree well 
with previously published all-sky galaxy counts, suggests 
that the phenomenon of a rising luminosity density with 
increasing distance from us may be a general trend in 
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